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Abstract 

We study the zero-temperature behavior of the infinite-ranged Ising spin glass 
in a transverse field. Using spin summation techniques and Monte Carlo meth- 
ods we characterize the zero-temperature quantum transition. Our results are 
well compatible with a value v = \ for the correlation length exponent, z = 4 
for the dynamical exponent and an algebraic decay t -1 for the imaginary-time 
correlation function. The zero-temperature relaxation of the energy in the 
presence of the transverse field shows that the system monotonically reaches 
the ground state energy due to tunneling processes and displays strong glassy 
effects. 
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I. INTRODUCTION 



The purpose of this work is to present some results concerning the zero-temperature crit- 
ical and relaxational behavior of the Ising spin glass in the presence of a transverse magnetic 
field. While classical spin-glasses have been extensively studied during the recent years, the 
role of the quantum effects in the low-temperature regime are not so well understood. In 
particular, a large amount of work has been devoted to the study of the one-dimensional 
case ]I] and the mean- field theory M. These two limiting cases seem to capture one 
of the most relevant features associated to the quantum fluctuations, i.e. the presence of 
tunneling effects at zero-temperature. The effect of the transverse field is to allow the sytem 
to jump over the free energy barriers even at zero temperature. In this work we will focus 
our attention in the study of the zero-temperature critical behavior and on the nature of the 
relaxational dynamics. We have considered the infinite-range model where some analytical 
results can be obtained. The infinite-ranged model has been studied in several works. In 
particular, the phase diagram of the model has been obtained using spin summation tech- 
niques while Miller and Huse || have obtained the imaginary-time correlation function 
at the zero temperature critical point using a theoretical analysis. On the other hand recent 
numerical work reveals that the Monte Carlo method can yield good estimates of the 
critical exponents associated to the quantum transition using finite-size scaling techniques. 

Our purpose is two fold. First we want to show how the Monte Carlo technique used in 
011 is a powerful tool in order to determine the critical point and the critical exponents in 
the mean-field case. This will be done comparing the results obtained using finite-size scaling 
and numerical spin summation methods. Once the critical field is obtained we will obtain 
the main critical exponents z and v and we will study the decay of the imaginary-time 
correlation function at the critical point. Unfortunately, our results are in disagreement 
with the theoretical prediction of Miller and Huse ||. Second, we will consider the role 
of the quantum fluctuations on the zero-temperature relaxational behavior of the model. 
While these last results concern the quantum infinite-ranged model we expect that our main 
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conclusions are valid also in the short-ranged case. 



II. THE MODEL 

The model we are interested is defined by the Hamiltonian, 

w = -£ v?°?- r Eof (i) 

i<j i 

where the {<jj; i = 1,N} are the Pauli spin matrices and T is the transverse field. The 
Jij are Gaussian distributed variables with zero mean and variance ■h. For T = the model 
reduces to the classical Sherrington-Kirkpatrick spin-glass model . It is well known |7j || 
that the ground state energy of the above Hamiltonian can be written as the free energy of 
a classical model with an extra imaginary-time dimension in the following way, 

EJT) = - lim lim M^lA (2) 

9K ' /3->ooM->oc Nf3 

where 

M N M 

Z eff = Tr ai exp(-(3H eff (T, M,{3))= £ exp(A ^ ^ J ij( j\ a) + 5£E o\a^ + C) 

<Ti=±l i<j t=l i=l t=l 

(3) 

and the spins Oi are classical variables which can take the values ±1. The parameters A and 
B are given by, 

M 

B= l -\og{coth{^)) (4) 

MN 1 . 2/3r 

C = \og{-sinh{——)) 

2 &v 2 v M it 

In the limit M — > oo the parameters A and S are highly anisotropic (the coefficient A 
goes to zero while B goes to infinity) . This makes extremely difficult to perform Monte Carlo 
numerical tests of the quantum model. It has been recently shown that it is better to 
work with a different Hamiltonian which nevertheless lies in the same universality class. To 



this end, we have considered the family of models with parameters A = (3a, B = (3 d , C = 0. 
Within this family of models the parameter (3 c i has the role of the inverse of a classical 
temperature (not to be confused with the real temperature) which controls the intensity of 
the quantum fluctuations. In some sense, this effective classical temperature -a- plays the 
role of a transverse field in the true model (|5J). In this way, the new Hamiltonian is a more 
isotropic one. Also, in case universality holds, we expect the mean-field critical exponents 
to be independent of the particular model considered. We have concentrated our attention 
in the previous models with n = 1 (model (a)), and n = 2 (model (b)) and we have studied 
them using the Monte Carlo method and spin summation techniques. While our Monte 
Carlo numerical results are consistent with the universality hypothesis we have discovered 
that model (a) is still hampered by strong Monte Carlo sampling problems while model (b) 
gives more confident results. 



III. SPIN SUMMATION RESULTS 

In order to apply the spin summation techniques we have analitically solved the previous 
model eq.(|]) using the replica trick with general coefficients A and B. The analytical solution 
of the infinite-range model has been already considered in the literature ||T4| , |T7i and here we 
will only remind the results. Applying the replica trick and performing the usual technical 
steps in the theory of spin glasses i.e. (introducing the order parameters and decoupling the 
different sites) one gets the effective free energy, 



M^)_ lim ^ = lim « (5) 



N(3 d n^o Nn(3 d *-*> n (3 d 

where (..) stands for average over the disorder and n is an integer which denotes the 
number of replicas. The saddle-point free energy A[Q, R] reads, 

A[Q,R] = ^(EY,(Q%) 2 + Y,Y,(K t ') 2 ) - log F[Q,R] (6) 

t,f a t£tf 

with 
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F[Q,R] =Eex P [BE^ +1 + j(EEW + EE*£V«Ol (7) 

The indices a,/3 = l,..,n stand for replica indices while the indices t,t' = 1,..,M run 
over the imaginary-time direction with periodic boundary conditions (i.e. a* m = al). The 
saddle point equations yield the order parameters Q and R, 

Q% = M) C = (M) (8) 

where the thermal averages (...) are done over the effective partition function defined in 
eq.(|7|). To solve the previous equation we impose the static condition (i.e. no dependence on 
the imaginary-time variables t, t') in the set of parameters Q while the R's are assumed to be 
no static but traslationally-time invariant, i.e. depend only on the difference of times t — t'. 
In order to determine the critical value of /3 c i it is enough to consider replica symmetry. In 
this case the order parameters assume the form = q, R l J' = R(t — t') an the free energy 
reads, 

f3f = ^-E( R ) - -s-(l"2?) - - / — - T e--log6(x) (9) 

4 t+t , 4 I J -co (2tt)2 

where the function Q(x) is given by, 
e(x) = E>xp(~(x,a))=5>xp(fi]>>V+ 1 + (A 2 g)^$>< + ^(R* - q)o*o*) (10) 

a* cr* t t t+t' 

and the order parameters q and Rit — t') can be obtained solving the equations, 
„_ f°° dx -^/ S^ t exp(S(a:,a)) x2 

^"^-i-oo^ol 6 2 1 e(i) ) (12) 

We have numerically solved the previous non linear equations for the models (a) and (b) 
at different values of M ranging from 2 to 15. Similarly as done in we have extrapolated 
the different parameters q and R(t — t') to the M — > oo limit. We have found that a second 
degree polynomial in i yields very stable and good results. In case of model (a) we found 
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a phase transition at T^' = 2.81 ± .01 while for model (b) we obtain T ci = 2.11 ± 0.01. 
The spin summation method yields the thermodynamic quantities with good precision but 
is inadequate to obtain the critical exponents at the transition. 



IV. MONTE CARLO RESULTS 

In order to characterize the quantum critical point we have done Monte Carlo (MC) 
numerical simulations of class of models (a) and (b). While model (a) displays strong Monte 
Carlo sampling problems (and needs a lot of computational time) the model (b) yields the 
critical behavior with modest computational effort. Note that model (a) corresponds the 
case considered in references In what follows, and otherwise stated, we will present 

numerical results for model (b). In order to simulate the system described by eq.(|3|) we 
consider M coupled systems along the time direction with the same realization of disorder. 
To increase the speed of the computations we have considered the case of discrete couplings 
Jij = ±— which yields the same behavior in the large N limit as in the case of a Gaussian 
distribution of couplings. We have simulated two different replicas {aj,T*;i = 1, ..,N;t = 
1, .., M} of the system eq.(|5|) with the same disorder realization. The main quantity we are 
interested in is the spin-spin overlap 

i N M 
i=l t=l 

which yields the spin-glass susceptibility, 



Xsg = N((q*) - (q) ) (14) 

Following J!]!! we consider the Binder parameter for different values of iV and M. This 
adimensional parameter measures the Gaussianity of the statistical fluctuations and is de- 
fined by, 



» = 5P-#>1 < 15 > 
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In the vicinity of the critical point the spin-glass susceptibility eq. (|14D and the Binder 



parameter eq.(|l5D are expected to scale with the size of the system N and the temporal 



dimension M in the following way, 

Xsg = N?x(N{T - T c y, N/M r ) (16) 

g = g{N{T-T c )\N/M r ) (17) 

where x, g are scaling functions and p, q, r are mean-field exponents related to the expo- 
nent v and the dynamical exponent z. [] 

Now we face the problem that the finite-size scaling depends on two variables N, M. 
As noted in [[| the phase transition is signalled by the behavior of the parameter g as a 
function of N and M. For large values of M the system behaves as a one-dimensional 
system and for small values of M the system behaves as the classical SK model. Then the 
Binder parameter ( |ToD is expected to go to zero for large and small values of of M. At 
intermediate values of Ma maximum for g is expected. Above the critical temperature the 
system becomes disordered and the value of g associated to that maximum decreases with 
N. Below T c it increases with N since the system tends to order. At the critical point T = T c 



the maximum value of g is constant with N. According to eq.(|T^) the scaling with N of the 
value of M corresponding to the position of maximum determines the mean-field exponent 
r. The previous criterium yields the critical temperature with very good precision. We find 
T = 2. 11 ±.01 in agreement with the results that we obtained in the last section. Our results 
for the spin-glass susceptibility \ and the Binder parameter g are shown in figures 1 and 2 
at T = 2.11. The values of N we studied cover the range N = 32 — 160 with 5000 samples 
in each case. We have observed that small values of N (in fact , less than iV ~ 50) are 
affected by strong subdominant corrections to the critical behavior. The reason is easy to 



1 This exponent z should not to be confused with the dynamical exponent asociated with the 
critical-time dynamics in classical systems. 
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understand since in the model we are studying the maximum of g is located at quite small 
values of M (for instance, at N = 32 the position value of M where the g has its maximum 
is located at a value less than 2 which is certainly very small). 

Larger values of N (we only show data for iV larger than 64) allow to extract the values 
of the critical exponents. The exponents p, q, r can be derived as a function of v and the 
dynamical exponent z. These are given by, v = y, zv = - which yield 7 = 2u. The 
numerical results for g show that the exponent r = | fits very well the scaling of the 
function g at the critical point. The fit of the spin glass susceptibility as a function of the 
temperature in the region of scaling M = 0.42iV5 is shown in the inset of figure 2 and is quite 
consistent with q — |,p = | which yields v = | and 7 = ~ as predicted within the Gaussian 
approximation [|I3| . Unfortunately it is difficult for us to conclude, from the numerical data, 
on the exact value of the exponent z. Our best fit reveals r = |, z = 3 which yields 
ft = |. But it is very plausible that these exponents are an artifact of the aforementioned 
subdominat finite M corrections. On the light of these considerations we expect that the 
canonical exponents r = | , z — 4 (which would also yield (3=1) are the correct ones. 
These are the values of the exponents used to scale data in figures 1 and 2. To definitely 
conclude on this point we should explore laregr sizes. But this is a very difficult task due 
to the long-ranged nature of the model we are studying which makes simulations very much 
time consuming. It is interesting to note that the critical value of g ( g c = Max(g(N, M, T c )) 
) is close to 0.056 and smaller than the values obtained in two and three dimensions 0|| 
as expected. As previously said we have also performed numerical simulations of model (a) 
which shows a critical value of g of order 0.07 slightly higher than that of model (b). But in 
this case we have not been able to make the data for g to collapse in a single universal curve. 
We are suspicious that strong Monte Carlo sampling problems are the reason for such bad 
results. This is presumably related to the value of B in the critical point which is higher in 
model (a) than in model (b). This implies stronger anisotropic interactions in the first case. 

Recently Miller and Huse have obtained the imaginary-time correlation function at the 
critical point using a theoretical analysis . Our mean-field exponents are in disagreement 
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with their results. At the critical point they obtain, 

C(t) = ~ t~ a (18) 

with the value a = 2. Figure 3 is a check of the theoretical prediction by Miller and 
Huse for the imaginary-time correlation function at the critical temperature T = 2.11 . 
Simulations have been done for a large system N = 2272, M = 20 such that it is in the scaling 
region where we expect the g(N, M, T c ) takes its maximum value. We have carefully checked 
that the system is in thermal equilibrium and data has been averaged over 8 samples. The 
results for the decay of the correlation function eq.(0) yields an exponent a ~ 1.2 consistent 
with the exponent a = which ranges from 1 to 7/6 depending if z = 3 or z = 4. Note 
that the decay of the imaginary-time correlation function eq.(^) is quite sensitive to how 
much close we are to the critical region. Obviously, if we are not precisely in the critical 
region we expect the system to be slightly more disordered and the correlation function to 
decay faster. In any case the fitted value 1.2 is an upper limit to the true exponent a which 
we find natural to be 1 and then z = 4. It is not clear to us how the predicted exponent 
a = 2 can fit the numerical data. 



V. ZERO-TEMPERATURE RELAXATION AL DYNAMICS 

Once we have characterized the zero temperature quantum transition we want to present 
some results concerning the real-time dynamical behavior of the quantum model at zero 
temperature. We face the problem of defining a reasonably real-time dynamics for a quantum 
system. We have considered the simple possibility that real time Monte Carlo dynamics is 
an appropiate tool to explore the slow dynamic process in the presence of tunneling effects. 
In the classical case (zero transverse field) we already know that the relaxation at zero 



temperature of the system stops whenever it founds a metastable state [Hi]. Because the 
dynamics is non ergodic in the classical case (the system cannot jump over energy barriers) 
then the system cannot reach the ground state energy. When a transverse field is applied 
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the system can jump over energy barriers allowing for a new type of relaxation. In order 
to investigate this point we have considered the relaxational dynamics of the true quantum 
model of eq.(|3]) with the coefficients A, B, C given in eq.(^|) at very low temperatures as a 
function of the transverse field. Concretely we are interested in the behavior of the model 
for large (3 in the limit M — > oo with jj as much small as possible 

In this limit the Hamiltonian eq.(f|) is strongly anisotropic, the coeflicent A goes like 
while B is much larger and goes like log(M). The total energy in eq.@ can be decomposed 
in two parts plus a configuration independent constant C: E = AEj + BEp + C where 
Ej is the sum of all interaction energies in the different imaginary-time slices and Ep is a 
nearest-neighbour ferromagnetic interaction between spins in the different imaginary-time 
slices. Our main quantity of interest is the relaxational behavior of the interaction part Ej 
as a function of time. We will show that the dynamical evolution of the system is the same 
if the Monte Carlo time is rescaled by the factor (^) 2 - This is a natural result since the 
parameters A and B of the effective Hamiltonian of eq.(|7|) are only a function of that ratio. 
Note that in the limit M — > oo the relaxation of the energy Ej is extremely slow with time 
(because the main contribution to the full energy in the Hamiltonian eq.(£J) is due to the 
ferromagnetic term Ep). This clarifies the apropriate regime of parameters (3 and M in 
which the zero temperature relaxation of the model is defined. Moreover, depending on the 
values of (3 and M one is considering, it unambigously determines a diferent region of the 
real dynamical time which is explored. 

We performed two kind of experiments. We have studied zero temperature dynami- 
cal relaxations at a fixed transverse field. We have considered the model at different low 
temperatures and different values of M such that yield nearly the same thermodynamic 
properties. In figure 4 we show the relaxation of the energy Ej as a function of Monte 
Carlo time for different values of M and (3 such that the ratio 4f is small. The simu- 



2 Note that in eq.@ the limit (3 — ► oo is performed after the limit M — > oo 
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lations were performed for two different sizes N = 320, 640 finding the same qualitative 
results. We studied several different ratios jj| ranging from 0.2 to 0.001. The explored tem- 
peratures were T = 0.1,0.02,0.01, deep in the zero-temperature region, and the values of 
M = 50, 100, 500, 1000. Relaxations were studied with a small transverse field T = 0.1 (the 
critical values of the transverse field is close to 1.6 0). In order to make the relaxation curves 
collapse in a single curve we have rescaled the time by the factor (^) 2 . It is interesting to 
observe that the energy Ej decreases to a value close to —0.76 which is the expected value 
in the classical SK model at zero temperature at first order of replica symmetry breaking 



| 15| . Note also that the energy Ej decreases with time but it can fluctuate and increase 
due to the tunneling effects in the presence of the transverse field. Indeed we have clearly 
appreciated this effect especially in the large-time regime. 

Another interesting aspect of the quantum model we are considering concerns its glassy 
properties due to tunneling effects. The transverse field controls the intensity of quantum 
fluctuations and we expect strong hysteresis effects as the transverse field is varied. This is 
shown in figure 6 where we plot the relaxation of the energy Ej at three different cooling- 
heating rates as a function of the transverse field T []. The cooling rate is defined by the 
number of Monte Carlo steps per temperature step (AT=0.05 in figure 6). Hysteresis curves 
for different values of M and (3 collapse in the same curve once the cooling rate is appropiately 
scaled by the time factor (jj) 2 - The area enclosed in the hysteresis curves decreases as the 
cooling-heating rate decreases very similar to what happens in real glasses. 

VI. CONCLUSIONS 

In this work we have studied the zero temperature behavior of the infinite-range quantum 
Ising spin glass in a transverse field. In particular we have studied the critical properties at 



3 In our case the parameter which is varied is the transverse field and not the temperature as in 
real glasses 
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the quantum transition point and the relaxational behavior as a function of the transverse 
field. 

Concerning the static properties we have studied an effective model (the so called model 
(b)) which is expected to be in the same universality class as the original quantum model 
eq. (|3],|5]) . Also this effective model does not present strong Monte Carlo sampling problems 
and gives enough confident results. Even though our results for model (b) show strong 
finite M corrections for small sizes, our data is in agreement with the mean-field quantum 



exponents v = 1/4, (3 = 1,7 = 1/2, z = 4 JT3|. Unfortunately we have not been able 



to corroborate the prediction of Miller and Huse for the imaginary time autocorrelation 
function eq.fllBD where a = 2. This is the result expected for a dynamical quantum exponent 
z = 2 which we definitely rule out from the analysis of the data shown in figures 1 and 3. 
In particular, numerical data shown in figure 3 reveals an exponent of a ~ 1.2 which should 
be a little bit lower if we are not precisely within the scaling region. The value a = 1 seems 
us the natural exponent compatible with our numerical results. This is an interesting point 
which deserves further investigation. Unfortunately it is very difficult to go to larger sizes 
since we would need much more computing time. 

Concerning the dynamical properties of the model we have investigated the zero temper- 
ature relaxational behavior of the model. We have found that the quantum model eq.([|[5|) 
in the zero temperature limit (3 — > 00, with — > shows a universal relaxational behav- 
ior when the Monte Carlo time is rescaled by the factor (^) 2 . For a low transverse field 
we have observed that the universal curve for the interaction energy Ej monotonically de- 
creases approximately to the static value predicted in the classical SK model at first order 
of replica symmetry breaking (obviously there are small corrections due to the finite value 
of the transverse field). Because the effective model (|3|||) mainly depends on the ratio jj| we 
expect that similar conclusions about the dynamical behavior of the infinite-ranged model 
are also valid in the short-ranged case. We have also observed the glassy features of the 
model by studying the hysteresis effects as a function of the cooling-heating rate variation 
of the transverse field. The results shown in figure 5 indicate a dynamical behavior of the 
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model reminiscent of that observed in real glasses. In the presence of a transverse field the 
system can jump over energy barriers due to tunneling effects. Then, at zero temperature, 
the system is not constrained to remain forever in a metastable state. It can be instructive 
to speculate if this jumping of the system over the energy barriers corresponds to some 
kind of activated processes in classical glassy models. This interesting point deserves further 
investigation. 
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Figure Captions 

Fig. 1 Binder parameter g(N,M) in model (b) at T = T c = 2.11 for different sizes N = 
64, 96, 128, 160 as a function of M/iV§. 

Fig. 2 Spin-glass susceptibility x{N, M)/N% in model (b) at T = T c = 2.11 for different sizes 
N = 64, 96, 128, 160 as a function of M/7V§. The inset shows the x( N , M)/N* as a 
function of the temperature for N = 32, 96, 160 for values of N, M where the g takes 
its maximum value. 

Fig. 3 imaginary-time correlation function in model (b) for N = 2272, M = 20 in the scaling 
region averaged over 8 samples. The fit is of the form C(t) = A/t a + A/ (20 — t) a with 
the best fit parameters a = 1.2, A = 0.3. 

Fig. 4 Relaxation of the energy Ej with T = 0.1 for different ratios (3/M as a function of the 
rescaled Monte Carlo time t* = t(jj) 2 . 

Fig. 5 Hysteresis cycles of the energy Ej at three different cooling rates r (dotted line r = 100, 
dashed line r = 10, continuous line r = 1). 
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